Astronomy & Astrophysics manuscript no. 11839 


June 17, 2009 


(DOl: will be inserted by hand later) 





o 
o 

(N 

1—5 



Oh; 
O 

(N 
> 
oo 
m 

O 
On 
O 

> 



Interactions of UHE cosmic ray nuclei with radiation during 
acceleration: consequences for the spectrum and composition 



D. Allard^ and R.J. Protheroe^ 

' Laboratoire Astroparticule et Cosmologie (APC), Universite Paris 7/CNRS, 10 rue A. Domon et L. Duquet, 75205 Paris 
Cedex 13, France. 

^ Department of Physics, The University of Adelaide, SA 5005, Australia. 

Received date; accepted date 

Abstract. 

In this paper, we study the diffusive shock acceleration of cosmic-ray protons and nuclei, taking all the relevant interaction 
processes with photon backgrounds into account. We investigate how the competition between protons and nuclei is modified by 
such acceleration parameters as the acceleration rate, its rigidity dependence, the photon density, and the confinement capability 
of the sources. We find that protons are likely to be accelerated to higher energies than nuclei in the case of interaction-limited 
acceleration processes, whereas nuclei are accelerated to higher energies than protons for confinement-limited acceleration. 
Finally, we discuss our results in the context of possible astrophysical accelerators and in the light of recent cosmic-ray data. 
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1. Introduction 

[ By plotting magnetic field vs. size of various astrophysical ob- 
' jects, Hillas d 19841 1 identified possible sites of acceleration of 
UHE CR based on whether or not the putative source could 
contain the gyro-radius of the accelerated particles, and on 
the likely velocity of scattering centers in these sites, and nar- 
rowed the field of possible sources to radio galaxy lobes and 
^ galaxy clusters. This issue has been re-examined recently in 
■ the context of diffusive shock acceleration (DSA) (Protheroe 
' 12004 ') taking account of energy losses (synchrotron) and in- 
teractions (Bethe-Heitler and pion photoproduction) which can 
cut the spectrum off, and so provide additional constraints. It 
was found that with the maximum possible acceleration rate 
for a given magnetic field the spectrum is cut off by pion photo- 
production if the magnetic field is below ~3 juG, and for higher 
fields the spectrum is cut off by synchrotron losses. 

However, for an acceleration rate 10 times the maximum 
possible, pion photoproduction cuts off the spectrum for fields 
below ~100 juG and synchrotron losses cut off the spectrum 
for higher fields. In order to accelerate protons above lO^^eV 
the ideal sources were found to have acceleration region sizes 
(magnetic field) between 1 kpc (~1 mG) and 1 Mpc (~1 fiG) 
such that the gyro-radius can be contained within the accel- 
erator These parameters allow the spectrum to extent beyond 
10^'^eV where it will be cut off either by synchrotron losses (for 
the higher magnetic fields in this range) or by pion photopro- 
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duction (for the lower fields in this range). This makes radio 
galaxies and clusters of galaxies the probable sites of UHE CR 
acceleration (see Fig. 6 of Protheroe 120041) in agreement with 
the conclusion of Hillas. However, if the whole acceleration 
region is moving relativistically, Doppler boosting in energy 
of neutrons (produced by pion photoproduction) escaping from 
such sources allows blazar jets (Protheroe et al. 120031 1 and GRB 
( Waxman 120061 and references therein) to remain possible sites 
of UHE CR acceleration, and in these cases proton synchrotron 
losses are likely to cut off the spectrum. Large-scale cosmic 
shocks in the Universe may also be sites of UHE CR accelera- 
tion (Kang & Jones 2002 Inoue et al. 120071) as well as colos- 
sal magneto-hydrodynamic fossil AGN structures (Benford & 
Protheroe 120081 ). AGN jets and cocoons surrounding AGN jets 
(Schopper et al. 120021 Norman et al. '1995 ) and hot spots in 
radio lobes of powerful radio galaxies (Rachen & Biermann 
119931 ). In these cases pion photoproduction, the size of the ac- 
celerator or the time available for acceleration may cause the 
cut-off in the spectrum of protons at acceleration. In a recent 
review, Biermann et al. (2()08) concluded that the most promis- 
ing contenders are radio galaxies and GRB. 

Apart from the nature of the cosmic accelerators, the com- 
position of the accelerated particles at the highest energies is a 
key parameter in understanding the origin of the extragalactic 
cosmic-rays and the physical conditions at play at the sources. 
To determine the spectrum and composition of accelerated par- 
ticles at an astrophysical source the injected composition at the 
shock and the interactions within the acceleration region should 
be considered. Satellite measurements of low energy galactic 
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cosmic -rays show that the cosmic -ray composition is deficient 
in H with respect to interstellar matter, and the relative abun- 
dances of the various species are believed to be related to their 
volatility (Meyer et al., 1997). It is of course possible that such 
enrichment takes place also in extragalactic sources. 

In the present paper, we take into account all the relevant 
interactions of protons and nuclei with ambient photon back- 
grounds and calculate the spectrum and composition resulting 
from diffusive shock acceleration for various compositions in- 
jected at the shock. Using a box model of DSA (Protheroe 
I2004I I. the SOPHIA event generator for pion photoproduction 
(Miicke et aI.. I2000l l and the photon-nucleus interaction event 
generator developed in Allard et al. (2005 2008), we study the 
dependence of the composition and spectrum of accelerated 
cosmic rays on the acceleration rate, the rigidity dependence 
of the diffusion coefficient, the confinement within the accel- 
eration site, and the density of the photon backgrounds. The 
paper is organized as follows, in the next section we introduce 
the relevant interaction processes of protons and nuclei with 
photon backgrounds. In Sect. 3 we describe the box model for 
DSA we use in our calculations. In Sect. 4 we present the spec- 
tra and composition resulting from the acceleration processes 
under different hypotheses about the injected composition and 
physical parameters for the accelerators. Finally, in Sect. 5 we 
discuss our results. 

2. Interactions of nucleons and nuclei with 
radiation and their modeling 

Protons and nuclei propagating through the intergalactic 
medium interact mainly with cosmic microwave background 
(CMB), infrared, optical and ultraviolet (IR/Opt/UV) photon 
backgrounds. These interactions produce features in the prop- 
agated UHE CR spectrum such as the "GZK cutoff" (Greisen, 
11966 ; Zatsepin & Kuzmin. 11966 ). Pions produced in the same 
interactions generate the cosmogenic neutrino flux (Berezinsky 
& Zatsepin, 1 19691 ). During the acceleration process, these inter- 
actions are also expected to take place at a rate depending on 
the ambient photon density and can also produce spectral and 
composition features in the cosmic-rays as well as secondary 
particles (neutrinos and photons), as discussed by Protheroe 
(120041 ) in detail for the case of proton primaries. 

In the following, we will use the recent estimate of Stecker 
et al. (120051 ) for the intergalactic background to model the pho- 
ton backgrounds in the infra-red, optical and ultra-violet win- 
dow. We re-scale its density with different scaling factors to ac- 
count for the higher density IR/opt/UV backgrounds that could 
be expected in astrophysical sources (see Fig.[T]). 

Bethe-Heitler pair production interactions with the various 
photon backgrounds are treated as a continuous energy loss 
process. We use the cross sections, inelasticities and the mass 
and charge scaling of the attenuation length for nuclei provided 
by Rachen (119961 ). If the energy of the background photons ex- 
ceeds ~145 MeV in the nucleon rest frame, protons and neu- 
trons can interact through the pion photoproduction process. 
This high inelasticity process is better treated by Monte Carlo 
calculation, and our Monte Carlo code uses the SOPHIA event 
generator (Miicke et al., 2000) allowing us to treat accurately 
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Fig. 1. Radiation field of the CMB compared with the 
IR/Opt/UV from Stecker et al. (12005b with various scaling fac- 
tors. Also shown, as will be discussed in Sect. 4, is the radiation 
field expected at 200 kpc along the jet from the flat-spectrum 
radio quasar (FSRQ) 0208-5 12 based on its spectral energy dis- 
tribution given in fig. 5 of Tavecchio et al. (120021 ) . 



the various interaction channels (direct pion production, reso- 
nances, multi-pion production) and their branching ratios (see 
Rachen [1 9961 ). The energy dependence of the mean free path 
of protons for photopion production with CMB and IR/opt/UV 
photons is shown in Fig.|2^. Also shown is the mean free path 
for a IR/opt/UV background 100 times denser. 

The interactions experienced by nuclei with photon back- 
grounds differ from those of protons. In addition to the pair 
production losses that result in a decrease of the Lorentz factor 
and the rigidity of the nucleus, one must consider the photo- 
disintegration (also called photo-erosion) processes that leads 
to the ejection of one or several nucleons from the nucleus. 
Different photo-erosion processes dominate the total interac- 
tion cross section at different energies (Puget et al., 119761 ). 
The lowest energy disintegration process is the Giant Dipole 
Resonance (GDR) which results in the emission of one or two 
nucleons and a particles. The GDR process is the most rele- 
vant as it has the highest cross section and the lowest thresh- 
olds, between 10 and 20 MeV, for all nuclei. For nuclei with 
mass A > 9, we use the theoretical GDR cross sections pre- 
sented by Khan et al. (120051) . which take into account all the 
individual reaction channels (n, p, 2n, 2p, np, a, . . . ) and are 
in better agreement with data than previous treatments. For nu- 
clei with A < 9, we use the phenomenological fits to the data 
provided by Rachen (1996). Between ~30 MeV in the nucleus 
rest frame and the photopion production threshold the quasi- 
deuteron (QD) process becomes comparable to the GDR, and 
dominates the total cross section at higher energies. 

Pion photoproduction (baryonic resonance (BR) process) 
by nuclei becomes relevant above 150 MeV in the nucleus rest 
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Fig. 2. (a) Left: Energy dependence of the mean free path for pion photoproduction of cosmic -ray protons in the intergalactic 
background radiation fields, and for the case where the IR/Opt/UV background is 100 times higher, (b) Right: Lorentz-factor 
dependence of the mean free path for photo-disintegration for four different nuclei (He, O, Si, Fe) in the same photon backgrounds 
as in the left panel. 



frame, i.e., ~ 5 x 10^' eV in the lab frame for iron nuclei in- 
teracting with the CMB. We use the parameterization given by 
Rachen (' 19961 ) in which the cross section in this energy range is 
proportional to the mass of the nucleus. The basis for this scal- 
ing that of the deuteron photoabsorption cross section which 
is known in great detail, although nuclear shadowing effects 
are expected to break this scaling above 1 GeV. It is impor- 
tant to note that pion photoproduction cross sections for nuclei 
are different from the free nucleon case. In particular, in nu- 
clei the baryonic resonances heavier than the first A resonance 
are far less pronounced than for nucleons, and the cross sec- 
tions can not be simply derived from the free nucleon case. We 
also follow Rachen (119961 ) for the treatment of secondary par- 
ticle production, nucleon multiplicities, energy losses, and the 
probability of absorption of the produced pion by the parent 
nucleus. 

The Lorentz-factor dependence of the total mean free path 
(all photo-disintegration processes and all channels) for photo- 
disintegration for four different abundant species is shown in 
Fig. I2J3 for the same photon backgrounds as Fig. [2^. One can 
see that the interaction thresholds occur at similar Lorentz fac- 
tors for all the species, and that the mean free paths scale ap- 
proximately linearly with the mass as a consequence of the 
mass scaling of the cross sections. Photo-disintegration is dom- 
inated by the GDR process over the whole energy range, except 
for the second minimum at the highest energies which is due to 
the BR process (see Fig. 3b of Allard et al. (120061) for the con- 
tributions of the various processes to the mean free path). 

3. Diffusive shock acceleration 

Diffusion of energetic particles in magnetic fields depends on 
their magnetic rigidity p = pcjZe, where p is the the parti- 



cle's momentum and Ze its charge. In DSA the rigidity gain 
and escape rates depend on the diffusion coefficient which is 
usually assumed to have a power-law dependence on rigidity. 
Kip) oc p*, with the exponent depending on the nature of the 
turbulence present in the magnetic field: 5=1/3 (Kolmogorov 
spectrum), 1/2 (Kraichnan spectrum) or 1 (completely disor- 
dered field). The diffusion coefficient is lowest, and the accel- 
eration fastest, for the case where the magnetic field B is highly 
disordered, and it can approach the Bohm diffusion coefficient. 
Kg - \^RgC where Rg - p/cB is the gyro-radius. In this case 
Jokipii (19871) showed that if the diffusion coefficient parallel 
to the magnetic field direction is a factor 77 > 1 times the Bohm 
limit, i.e. 

1 



(1) 



then the diffusion coefficient perpendicular to the magnetic 
field is expected to be approximately 

kii 



\+Tf- 



(2) 



provided that 77 is not too large (values in the range up to 10 
appear appropriate). 

Under this assumption, and for the case of DSA at a strong 
plane parallel shock (i.e. with magnetic field B parallel to the 
shock normal) with shock speed u\, this leads to a rigidity gain 
rate 

1 dp 

rgmnip) = 7 — 



p dt 



3u\B 



(3) 



and for a perpendicular shock (i.e. with magnetic field perpen- 
dicular to the shock normal) the rigidity gain rate would be a 
factor ~ 2.577(1 + 77-) higher (Protheroe 1998] l. 

In general, for particle acceleration by electric fields in- 
duced by the motion of magnetic fields B (including those at 
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astrophysical shocks), the maximum rigidity gain rate of rela- 
tivistic particles can be written 



(4) 



where ^(p)< 1 is the acceleration rate parameter which depends 
on the details of the acceleration mechanism, with ^=1 cor- 
responding to the highest posible gain rate for any mecha- 
nism. For DSA ^(p) will depend on the shock velocity, mag- 
netic field alignment and diffusion coefficients. For example, 
for the parallel shock case with mi = 0.1c and 77 = 10 one finds 
^(p) < 1.5 X 10 and for the equivalent perpendicular shock 
case ^(p) < 0.04 (see Protheroe ll998l) . Detailed and rigorous 
treatments of DSA are given in several review articles (Drury 
[1983 Blandford & Eichleril987. Berezhko & Ki-ymsky 09881; 
see particularly the review by Jones & Ellison ( 119911 ), on the 
plasma physics of shock acceleration, which also includes a 
brief historical review and refers to early work. 

Protheroe and Stanev (119991) proposed a box model of DSA 
(based on that developed by Szabo & Protheroe 119941 ) which 
was able to reproduce the essential features of DSA in a sim- 
ple leaky-box scenario. Particles were injected into the "box", 
and while inside the box their energy increased at the aver- 
age rate that would apply given the physical parameters of the 
shock acceleration being modelled, i.e. the upstream and down- 
stream diffusion coefficients, the shock velocity and compres- 
sion ratio. Particles being accelerated were made to leak out 
of the box, again at the average rate that would apply given 
the physical parameters, and these particles represented the ac- 
celerated particles escaping downstream in the standard shock 
acceleration picture. Protheroe ( 12004!) incorporated improve- 
ments suggested by Drury et al. ( 119991 ) into the model, and 
used it to investigate cut-offs and pile-ups in spectra of accel- 
erated particles in the presence of energy losses or interactions, 
including the case of acceleration and propagation of UHE cos- 
mic ray protons in the presence of the CMB. It was found that 
although almost all information about the acceleration condi- 
tions was lost after propagation, as a result of the GZK cut off, 
significant information remained imprinted in the spectrum of 
cosmogenic neutrinos produced in the "GZK interactions". In 
the present paper we use the calculation scheme of Protheroe 
( I2004i ) as it is very convenient to implement in a numerical or 
Monte Carlo program to investigate the spectrum and compo- 
sition of UHE cosmic rays accelerated in the presence of the 
CMB and IR/Opt/UV radiation fields. 

As we have seen, the normalization of the gain rate depends 
in a non-trivial way on the magnetic field, turbulence spectrum, 
shock configuration and shock velocity. Different classes of 
source, and even different sources of the same class, will ac- 
celerate particles up to different maximum energies. We have 
therefore arbitrarily chosen to normalize the gain rate at a given 
cut-off rigidity to that necessary for protons to reach that given 
cut-off rigidity in a steady state situation in the presence of the 
CMB and IR/Opt/UV radiation. We have done this irrespective 
of the rigidity dependence of rgain(p) for two reasons: firstly be- 
cause we are interested in the maximum energies achieveable 
and in interactions occurring at or near cut-off, and secondly 
because it is only the gain rate close to the cut-off rigidity, not 



the rigidity-dependence of the gain rate, that determines the 
value of the cut-off rigidity, at least for a steady-state situation. 



The nominal cut-off rigidity for protons pfj^ occurs where 
the gain rate equals the loss rate in the intergalactic IR/Opt/UV 
background. 



(Pcu!) = '"io.ss(p; 



1). 



(5) 



'gainVPcut 

In Fig. [3] we show the energy dependence of the loss (or at- 
tenuation) length for cosmic -ray protons, xiofis(p,Z = 1), in 
the CMB plus intergalactic IR/Opt/UV background (continu- 
ous line) and in the CMB plus an intergalactic IR/Opt/UV back- 
ground 100 times higher (dashed line) (assuming the same pho- 
ton backgrounds as in Fig.|2]). We have added to Fig.[3]curves 
representing Xg^i^ip) = c/rgai„(p) for different values of p^^J 
(assuming 6=1/3 and 6=1). Using the value of p^^'l in the CMB 
plus intergalactic IR/Opt/UV background to define the normal- 
ization of the rigidity gain rate we have 



'gain 



(P) 



Pcut/ 



'"loss(Pcm, 



1) 



Pcut / 



(6) 



We have chosen for p*^,' values in the range 10'*^ V to 
10^^ V, and we now discuss for what magnetic fields and ac- 
celeration rate parameters it is feasible for protons to achieve 
these cut-off rigidities in the CMB plus IR/Opt/UV radiation. 
For p[''J, = 10'^ lO^o, 10-2 and 10-^ V, the corresponding 
values of ^(p<^,5)B are 1.5 x 10-'^ 6.9 x lO"''* , 6.9 x lO^'^ 
and 1.6 x 10"'" T. However, for a given value of the product 
^(Pajt^B, not all vales of B are possible as protons will lose en- 
ergy also by synchrotron radiation, and this actually becomes 
more important than pion photoproduction losses at high rigid- 
ity. If the proton spectrum is cut-off by synchrotron loss, then 




19 20 



22 



Fig. 3. Energy dependence of the loss length for cosmic -ray 
protons in the CMB plus intergalactic IR/Opt/UV background 
(continuous line) and assuming a IR/Opt/UV background 100 
times higher (dashed line). The loss length is compared to the 
gain length assuming three different values forp^^/, (10'^ lO^" 
and 10^2 V) and 6=1/3 (thin short dashed hnes) and 6=1 
(thin long dashed lines). 
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Fig. 4. Cut-off rigidity of protons accelerated in the presence 
of CMB radiation vs. magnetic field for ^-1, 0.04 and 1.5 x 
10 chain curve - cutoff due to Bethe-Heitler pair production 
and pion photoproduction; dashed curve - cutoff due to syn- 
chrotron; sohd curve - total. 

the cut-off rigidity is given by 

p,^J = 2xlO'V^'B"''' V (7) 

where B is in tesla. To illustrate this, we plot in Fig.|4]the cut-off 
rigidity vs. magnetic field taking into account of synchrotron, 
Bethe-Heitler pair production and pion photoproduction losses 
for three values of acceleration rate parameter For example, for 
the parallel shock case with M] = 0.1c and 77 = 10 for which 
^(p) a: 1.5 X 10""* it is possible for protons to reach rigidities 
above p - 10''^ V for magnetic fields in the range 5 x 10"'^ to 
5 X 10"* T, and to reach rigidities above p - 10^" V fields in 
the range 10""^ to 5 x 10"*^ T are required with the maximum 
rigidity possible for that case being 1.5 x 10^" V achievable for 
10"^ T. For the equivalent perpendicular shock case for which 
^(p) «; 0.04, p = 10^' V can be reached for magnetic fields in 
the range 10"'" to 2 x 10"^ T, with the maximum rigidity pos- 
sible being 5 x 10^' V achievable for 10"^ T. Higher rigidities 
require relativistic shocks or different mechanisms, but proton 
rigidities above 4 x 10^^ V are are not possible without Doppler 
boosting of neutrons produced in pion photoproduction within 
a relativistically moving source such as an AGN jet or GRB. 

In Fig. |5] the energy dependence of the loss length for He 
(Fig.|5t) and Fe (Fig.|5}3) is plotted (assuming the same photon 
backgrounds as in Fig.|2|i and is compared to the energy depen- 
dence of gain lengths for the same rigidity dependences as the 
curves in Fig. [3] As the abscissa refers to the energy, and the ac- 
celeration depends on rigidity, the curves for the gain length for 
a given p'^J are shifted in energy proportionally to the charge 
(i.e. by a factor of 2 for He and a factor of 26 for Fe). Although 
the loss length is not strictly relevant for nuclei because nuclei 
do not remain on the same energy loss curve when losing en- 
ergy, unlike protons, these figures show that the output of the 
acceleration process will crucially depend on the normalization 
of the background radiation and the gain rate and its rigidity 
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Fig. 5. As Fig.[3]but for (a) He (top), and (b) Fe nuclei (bottom). 

dependence, and that the "competition" between protons and 
nuclei will, at least for some combinations of parameters, be 
non-trivial. 

3.1. Box model approximation 

In the box model of shock acceleration, particles of magnetic 
rigidity po are injected into the acceleration zone, or "box", and 
while inside the box are accelerated at a rate r^otip) and escape 
from the box at a rate resc(p)- These two rates uniquely deter- 
mine the spectrum of accelerated particles, i.e. those escaping 
from the box. 

The average rate of change of rigidity of a nucleus of rigid- 
ity p and atomic number Z during shock acceleration in the 
presence of energy losing processes is 

^ = prtotip, Z) = p [rgain(p) - noss(p, Z)l (8) 

at total,Z ^ ^ 

where rioss(p,Z) includes both losses due to interactions such 
as Bethe-Heitler pair production that can be approximated as 
continuous energy losses, and pion photoproduction in which a 
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significant fraction of a proton's energy is lost per interaction. 
Hence, using Eq.|6]we have 



rtot(p,Z) = rioss(pJ.'^*,Z= 1) 



Pcut / 



- rioss (p,Z). 



(9) 



If there are no energy losses and no additional escape pro- 
cesses a power-law spectrum results, and the integral spectral 
index , (F-l), equals the ratio of the escape rate to the gain rate. 
In this case, the escape rate representing escape downstream is 
(F- l)?"gain(p), and for injection of A^o particles of rigidity po the 
expected differential rigidity spectrum of accelerated particles 
is 



cDo(p) = M)(r-i)p;;-y^. 



(10) 



The acceleration zone extends distances Li(p) = Ki(p)/ui 
and L2(p) = K2{p)/u2 upstream and downstream from the 
shock, respectively where ui and U2 are the upstream and 
downstream flow velocities in the shock frame. Adding an extra 
escape term r^™(P'Pmax) responsible for the "maximum rigid- 
ity" for example arising from the finite size of the accelerator, 
and an escape term rg||!°"'(p) resulting from energy losses caus- 
ing particles to "fall out of the box" gives the total escape rate 



resc(p) = (r-i)nos.s(p^^j,z=i) 



p 

^Pcut / 



+ C(P,Pm.x) + rfsc™V,Z). 



(11) 



For the case of diffusion (with diffusion coefficient k oc p'') per- 
pendicular to the shock normal to the edge of the acceleration 
region being responsible for the maximum rigidity. 



(P,Pmax) = (F- 



max 



l)nos.s(Pcur 



Z=l) 



A^max 
PPcuxJ 



For the escape rate resulting from energy losses 



fallout 



(P,Z) 



1 dL2 dp 



L{p) dp dt 



= S ^2 noss (P,Z), 



(12) 



(13) 



where L-(Lx+ L2), and we adopt (2 « 0.5 (Protheroe |2004| l. 

Interactions of relativistic particles with matter or radia- 
tion generally result in the interacting particle losing energy, 
and this energy lost is either given to the struck particle or 
photon (elastic collision) or used in production of secondary 
particles (inelastic collision). If the mean interaction length 
(mean free path) is Amtip), for ultra-relativistic particles the 
rate of interaction is r-miip) - clA\nt{p). If the mean inelastic- 
ity, i.e. average fraction of energy lost per interaction, is a{p) 
then the effective loss rate for these non-continuous losses is 
'^loss O') - o;(p)'"mt(p)- Of course, if a(p) «: 1 then one can ap- 
proximate the interactions as continuous energy losses. 

In the Monte Carlo simulations Bethe-Heitler pair pro- 
duction was treated as a continuous loss processes, but 
hadronic collisions of protons and neutrons as well as photo- 
disintegration of nuclei were treated by the Monte Carlo 
method. Neutrons produced, e.g. in pj nn^ or during photo- 
disintegration of nuclei {Ny N' + in + jp, where and A^' 
are the parent and daughter nuclei, and ; and j are integers). 



may decay inside the box into protons which continue to be 
accelerated, or they may escape from the box to decay outside 
into cosmic rays depending on the neutron's Lorentz factor and 
the physical dimensions of the acceleration region. Since the 
simulations are for a given p^^J, p^ax and 6, and the size of the 
box depends on the (unknown) normalization of the diffusion 
coefficients, we have performed simulations separately for the 
two extreme cases. 

In the simulation with A^o particles injected, a particle is 
injected at time t-Q with rigidity p(0) - po and statistical 
weight w(0) - Wo = 1 /Nq. Its subsequent rigidity, p(f), and 
weight, w(0, are determined after successive time steps Af 
chosen to be much smaller than the smallest time-scale in 
the problem. In each time step, first the rigidity is changed, 
p{t + At) = p(f)[l + Af racc(p)], and then the probability of es- 
caping in time At is estimated as Pesc = {1 - exp[-Af resc(p)]). 
Then w(f)^'esc particles with energy E are binned in the his- 
togram of accelerated particles for the particle's species, and 
the particle's weight is changed to reflect the fraction not es- 
caping, W{t + At) = W(f)(l - Pesc)- 

Next, the probability of interacting in time At is estimated, 
Pmt,A,z - {1 ~ exp[-Af rint(p, A,Z)]j, and a random number 
is generated to determine whether or not an interaction takes 
place. If an interaction does take place, the energy of the target 
photon in the proton or nucleus rest frame (NRF), e', is sam- 
pled from 

p(e') oc n(e')o-(e') for fjj,^ < e' < 2ye^^^ 

where n(e)de is the LAB-frame number density of target pho- 
tons with energy in the range e to (e -H de), cr(e') is the relevant 
cross section in the NRF, e^'j^^. is the NRF threshold energy for 
the process considered, y is the Lorentz factor of the proton or 
nucleus, and fmax is maximum target photon energy in the LAB 
frame. For ultra-relativistic projectiles, the photon's direction 
in the NRF will be anti-parallel to the LAB-frame direction 
of the projectile. In the nucleon case (either a primary proton 
or a nucleon resulting from photo-disintegration), pion photo- 
production interactions are simulated in the the NRF using the 
SOPHIA event generator, and the 4-momenta of all the pro- 
duced particles are Lorentz transformed to the LAB frame. If 
the leading nucleon remains a proton, its rigidity (p, before the 
interaction) is changed to p = p,(l - a), where a is the inelas- 
ticity for that collision. However, if it is a neutron we consider 
two extreme possibilities as discussed earlier: (a) it decays to a 
proton inside the acceleration zone and will continue to be ac- 
celerated, or (b) it escapes from the acceleration zone, decaying 
outside and adding to the pool of cosmic rays which have es- 
caped from the accelerator. Escaping protons, neutrons, nuclei 
and neutrinos are binned in energy. 

In the case of interactions of nuclei, the number of emitted 
nucleons and a particles is determined according to branching 
ratios of the various interaction channels. Whenever the BR 
process is chosen, we check whether or not neutrinos are pro- 
duced and calculate their energy according to the A resonance 
kinematics (see Allard et al. l2006l for more details) 

Next we take account of the additional escape resulting 
from the rigidity decreasing during the interaction. This is 



D. Allard and R. J. Protheroe: Interactions of UHE cosmic ray nuclei with radiation during acceleration 



7 



more important for pion photoproduction losses than for photo- 
disintegration in which the rigidity does not necessarily de- 
crease. The probability of the particle immediately escaping 
downstream due to falling out of the box is estimated. 



Prob. (escape, p; p) - 



Upi) 



1 



f2, (14) 



and a random number is generated to determine if this hap- 
pens. If the particle does escape, then w{t + At) particles with 
momenta p are binned in a histogram of accelerated particles, 
and a new particle is injected with rigidity po and weight wq. 
If the particle does not escape, then the particle's rigidity and 
weight are evolved through a new time step as described above. 
Finally, if the time since injection exceeds a maximum time 
which may represent the time during which the accelerator is 
active, factive, the acceleration ceases and a new particle is in- 
jected. In most cases we set to factive=10 Gyr, but we shall ex- 
plore shorter activity time scales. 

We note that the term rioss(Pcut'-^-l)' ^^^<i to normalize 
the rigidity dependence of the acceleration rate, represents the 
loss rate of a proton at a its cut-off rigidity p^^J in the extra- 
galactic photon background (i.e. the CMB and the inergalac- 
tic IR/opt/UV background). Using the parameter pJ,^J (i.e. the 
rigidity at which the loss rate would equal the acceleration rate 
for protons in an accelerator "bathed" in the intergalactic pho- 
ton background) as a reference will turn out to be very conve- 
nient for estimating the effect of different physical parameters 
(Pmax, ^ and the scaling factor of the IR/opt/UV background) 
on the spectrum of accelerated protons and nuclei. Our Monte 
Carlo results and the conclusions of our study will be presented 
in the following two sections. 



4. Results 

In this section, we describe the results of our calculations for 
different combinations of the above-mentioned acceleration pa- 
rameters. For each combination of 6, pmax and the scaling of 
the IR/opt/UV, we calculate the expected output spectrum for 
different values of p^.^* (10^°, 10^', 10^^ V) and various com- 
positions injected at the shock: (i) pure proton, (ii) pure iron, 
and (iii) a mixed composition similar to low energy Galactic 
cosmic rays (see Allard et al., 120051 and references therein). 
In Figs. l6l-[m we give the spectra of nuclei, showing contri- 
butions of different groups of elements, accelerated such that 
in the absence of interactions and losses the spectra would be 
proportional to p and in the present work we set the nominal 
spectral index to be T-l. The choice of other reasonable values 
of r would make only a minor difference to the relative com- 
positions of accelerated particles, and the relative contributions 
of secondary nucleons, photons, and neutrinos would be some- 
what higher (lower) for harder (softer) spectra, but this would 
not modify the main conclusions of the present work. In these 
figures we also show the flux of neutrinos plus antineutrinos 
(summed over all flavors) produced during acceleration, and 
for the case where neutrons are assumed not to decay inside the 
acceleration zone we give the contribution to the total escaping 
cosmic rays of neutrons produced during pion photoproduction 



and photo-disintegration of nuclei. Note that unless otherwise 
specified we assume that the source is active for 10 Gyr, and 
that the injection rate does not depend on time. 

Fig. |6] shows the accelerated cosmic -ray spectra as- 
suming Pmax - 10^^ V (corresponding to an infinite 
shock), Kolmogorov turbulence (diffusion coefficient rigidity- 
dependence index 6 = 1/3), a photon background equivalent 
to the extragalactic background, and neutrons decaying inside 
the acceleration zone. For this case, the acceleration process 
is interaction-limited because pn,ax is extremely high, and we 
see that there is no obvious scaling of the maximum energies 
of the various species with their mass or charge - this is con- 
sistent with our expectation based on Figs.[3]and|5] Indeed, al- 
though the interaction threshold for a given photon background 
will occur at energies energies more or less proportional to the 
mass number of nuclei, such a scaling is not observed for the 
various cases studied in Fig. |6l This is because the mean free 
paths also scale with the mass, which means that the energy at 
which the photo-disintegration mean free path becomes of the 
same order as the acceleration length scales in a non-trivial way 
with mass number A. However, heavy nuclei can reach higher 
energies than light and intermediate nuclei due to their higher 
interaction energy threshold and higher acceleration rate at a 
given energy. 

Once the acceleration rate becomes lower than the inter- 
action rate for a given component, a sharp cut-off takes place 
as the nucleus changes to a lower mass and lower energy. As 
can be seen in Figs. [3] and |5] the cut-off energy depends on 
p^yj and 6 (see below). Due to the very low values of the in- 
teraction length for nuclei with CMB photons, such cut-offs 
seem to be extremely difficult to avoid even for very rapid ac- 
celeration (corresponding to high values of pj:^!). The proton 
component also experiences a cut-off at energies depending on 
p*y'( because of the increase of the photopion interactions loss 
rate with energy and the subsequent possibility of protons to 
fall out of the box, the latter being more important for protons 
than for nuclei because nuclei mainly lose energy by ejecting 
nucleons. This cut-off is however smoother than in the case of 
nuclei due the different nature of the energy losses. In the case 
we study here, i.e. acceleration not limited by confinement and 
a source active during the Hubble time, protons can reach an 
energy slightly higher than heavy nuclei, their lower accelera- 
tion rate (at a given energy) being compensated for by a lower 
interaction rate. When the acceleration rate is high enough (e.g. 
p*yj > 10^' eV) most species of nuclei (even the lightest) can 
reach energies above 10^*' eV. 

For the case of a mixed composition or 100% iron there 
is a bump in the spectrum at the highest energies due to sec- 
ondary nucleons, and it starts quite abruptly at the energy at 
which the photodisintegration interaction rate becomes higher 
than the acceleration rate. Secondary protons and protons from 
decay of secondary neutrons in the box as assumed Fig|6]con- 
tinue to be accelerated and produce a tail in the energy spec- 
trum at the highest energies. 

The alternative case in which neutrons escape before decay- 
ing is equivalent to adding a supplementary escape term as ap- 
proximately half of the produced secondary nucleons leave the 
box without being further accelerated. In this case, protons also 
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Fig. 6. Spectra of the various nuclei emerging after acceleration together with the spectra of neutrons and neutrinos (sum of 
all flavors) escaping from the acceleration region. The initial composition is 100% H (left column), 100% Fe (middle column), 
mixed composition representative of that observed in Galactic cosmic rays (right column). The acceleration rate is such that 



p'^J = 10^° V (top row), 10^' V (middle row), and 10^^ V (bottom row). The IR/Opt/UV background is taken to be that at z=0, 

^max - 10^^ V. Other acceleration parameters are 



neutrons are assumed to decay inside the acceleration zone, K{p) oc p^^^ and p^a, 
as indicated. Curves are for atomic numbers as indicated in the legend, with the topmost solid curve being the total spectrum, 
middle solid curve being Z- 1 , and the bottom solid curve being for the total neutrino spectrum. 



have an additional probability of escaping after each pion pho- 
toproduction interaction equivalent to the probability of isospin 
flip (see Miicke et al., 2000). As a consequence, the cut-off in 
the proton component in Fig. |7]is sharper than for the case of 
neutron decay inside the box shown in Fig|6] In the 100% pro- 
ton case there is a bump in the spectrum due to neutrons which 
starts at the energy at which the pion photoproduction interac- 
tion rate becomes higher than the gain rate for a given value 
of p^^jj. When nuclei are injected there is also a second bump 
at lower energies corresponding to the interaction threshold for 
photo-disintegration. 



Case of a high ambient photon background 

As mentioned earlier, the photon background in astrophysi- 
cal sources could be significantly higher than in intergalactic 
space, and to mimic such stronger fields we have scaled the 
IR/opt/UV intergalactic background upward by various factors 
to investigate the influence of the photon density on the accel- 
erated spectrum and composition. We have kept the other pa- 
rameters the same as for Fig.|6] i.e. 5 - 1/3, pmax = 10^^ V and 
neutrons decaying inside the box, and show in Fig.[8]the result 
for an IR/opt/UV background 100 times higher than in inter- 
galactic space. As can be seen in Fig.[Tl such a background is 
comparable to the background expected in radio lobes located 
~200 kpc from a relatively powerful FSRQ and so can serve as 
a typical reference for what can be expected in some classes of 
source (see discussion below). 
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Fig. 7. As Fig. |6] except that neutrons are assumed to rapidly escape from the acceleration zone where they decay into protons. 
Curves are as in Fig.|6]except that the additional thick solid curve shows the contribution of neutrons. 



Figs. |3] and |5] show that such a high background will lower 
the energy at which the spectrum of a particular species cuts off 
as the balance between acceleration and energy losses/photo- 
disintegration losses occurs at a lower energy. We see in Fig. |8] 
that this effect is stronger for nuclei than protons, but that 
nevertheless most species can still be accelerated well above 
10'^ eV. We have found that this apphes even for IR/opt/UV 
backgrounds 1000 times higher than in the extragalactic space, 
but that when IR/opt/UV background is 10"* times higher than 
in extragalactic space protons and nuclei are not accelerated 
above lO'^ eV whatever the value of p[.^*. 



High photon backgrounds present at sources are usually 
thought to eliminate nuclei during the acceleration process giv- 
ing rise to a 100% proton composition. However, we find this 
not to be the case for any combination of the acceleration pa- 
rameters that we have explored. 



Influence of rigidity dependence of the diffusion coeffi- 
cient 

The energy dependence of the acceleration rate depends on the 
rigidity dependence of the diffusion coefficient k cc in the 
vicinity of the shock. As can be seen in Figs. [3] and |5] the 
power-law index 6 affects the competition between accelera- 
tion and energy losses. We show in Fig. |9] results for 6=1 
(Bohm scaling) while keeping all other parameters the same 
as used for Fig.[8]for which (5=1/3 (Kolmogorov turbulence) 
was used. We see that as result of the steeper dependence of the 
acceleration rate on energy the cut-offs in the spectra of the var- 
ious nuclear species are shifted to higher energies with the gap 
between protons and nuclei being reduced, and heavy nuclei 
being able to reach as high energies as protons even with the 
IR/Opt/UV background 100 times that in intergalactic space. 

We introduce here FiglTO]to illustrate some of the effects 
and dependencies that we have noted already, and to highlight 
them in a form which more clearly shows their origin. The ef- 
fect of changing the rigidity dependence of the acceleration rate 
noted above is even more pronounced if neutrons escape from 
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Fig. 8. As Fig.|6]except that the IR/Opt/UV is taken to be 100 times that at z-Q. 



the acceleration zone (see Fig. [TOf). The comparison between 
Figs.fTOb andflOb shows there to be a dramatic increase in the 
maximum energies of protons and nuclei in a very high radia- 
tion field (a factor lO** increase in IR/Opt/UV) when the rigid- 
ity dependence of the diffusion coefficient is changed from p'^^ 
(Kolmogorov) to p' (Bohm). 

Another consequence of a stronger dependence of the ac- 
celeration rate on p is a sharper cut-off of the proton component 
because the gain rate more rapidly decreases above p*^J and it 
is then very difficult to continue to accelerate protons once they 
start to interact, even if they do not escape from the source im- 
mediately. 

Confinement-limited acceleration 

So far we have only considered cases where the acceleration 
is limited by interactions which we model by setting pn,ax to 
a very high value, 10^^ V. However, astrophysical sources are 
of finite size and so are expected also to have limited capa- 
bihty of confining high rigidity cosmic rays, e.g. when their 
larmor radius becomes comparable to the size of the accelera- 



tion medium. Confinement-limited acceleration processes can 
be modeled by lowering the value of pmax- Fig. [TT] shows the 
result of our calculations assuming p^ax = 10'"^ V while keep- 
ing all the other parameters the same as in Fig. [8] (neutrons 
decay inside the box, 6=1/3 and a scaling factor of 100 for 
the IR/opt/UV background). In all cases, we see that the proton 
component cuts off around 10'^ eV due to the confinement lim- 
itation. With a rigidity cut-off one would expect that in absence 
of energy losses the cut-off energies would be proportional to 
the charge for the various species, but one can see in Fig.fTTIthat 
for this choice of the parameters this scaling only applies very 
crudely and only for the rapid acceleration case (p[.^* = 10^^ 
eV). The acceleration of nuclei still appears to be interaction- 
limited, but nevertheless heavy nuclei are always accelerated 
above lO^" eV when p<^J > lO^' eV and are accelerated to 
higher energies than the proton component. 

We find that for pmax = 10'^ V provided there is a steep 
energy dependence of the rigidity gain rate (6=1) the expected 
scaling of the cut-off energies with Z applies for the case of 
a photon background with the IR/Opt/UV 100 times that in 
intergalactic space (top panel of Fig. [TO]i and even with the 
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Fig. 9. As Fig. [8]except that K(p) oc p' . 



IR/Opt/UV 1000 times higher (Fig.[TOfe). The expected scal- 
ing also applies for lower photon backgrounds. 

For higher values of p,nax, e-g- 10^' V, there will be very 
little difference for nuclei to the p^ax - 10^^ V case (e.g. p|,^J = 
10^' V and 10^^ V in Fig.[8]l as interactions akeady cut off their 
spectra at rigidities below ZxlO^' eV. Only the proton spectrum 
will be confinement-limited reducing the gap between the cut- 
off energies of protons and nuclei. 

Influence of duration of source activity 

Some of the astrophysical sources that are candidates for the 
acceleration of cosmic-rays to the highest energies, e.g. FSRQ 
which according to unification schemes (e.g. Urry & Padovani 
d 19921 1 are the beamed counterparts of the powerful Fanaroff- 
Riley Class 11 (FR-11) radio galaxies, are transitory sources with 
a life time small compared to the age of the universe. The finite 
duration of the source activity (by activity we mean the period 
during which the source is able to accelerate cosmic-rays) does 
not necessarily impact on the spectrum and composition of ac- 
celerated cosmic rays. Indeed, if the acceleration time scale is 



small compared to the duration of source activity, there will 
essentially be no difference. An example is shown in Fig. [TOl . 
where the source parameters of Fig. [8] are used except for the 
source activity duration which is only 100 Myr We see that for 
a high acceleration rate (p|,^j = 10^^ V) the output is essen- 
tially identical to that shown in Fig. [8] with the nuclear com- 
ponents being almost indistinguishable between the two cases 
and the only major difference being that the high energy tail of 
the proton component is absent when the source has the shorter 
lifetime. 

The situation changes when the source activity duration is 
less than the time needed to accelerate protons at the highest 
energies. Fig.fTOh shows the spectrum and composition for the 
same parameters as Fig.[TOj except that the source activity du- 
ration is only 10 Myr. In this case, the proton component is even 
more affected due to its lower acceleration rate and the maxi- 
mum energies scale proportionally to the charge of the species. 
In fact the proton cut-off rigidity, at 6 x lO'^ V, is significantly 
lower than 10"" V which would naively be expected based on 
reading off p in Fig. |3]at x-3 Mpc (corresponding to 10 Myr) 
for the p[.yj - 10^^ V and x^ain oc p'/^ case. The reason for this 
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log(£' / eV) 

Fig. 10. Miscellaneous combinations of parameters: (a), (b) and (c) As Fig.|6]except that 5 - 1, Pmax - lO'^ V and only the mixed 
composition is shown, (d) Mixed composition, p^'^^J - 10^^ V, IR/Opt/UVxlOOOO, 6 - 1/3, pmax - 10^^ V and neutrons decay 
inside the box. (e) Same as (d) but with <J = 1. (f) Mixed composition, p[,^J = lO^^ V, IR/Opt/UVxlOOO, 6 = 1, p,™x = 10^^ V and 
neutrons escape from the box. (g) Same as (c) but with IR/Opt/UVxIOOOO. (h) Same as Fig.[8]for the mixed composition and 
p'yj = 10^^ V except that the source activity duration is 10 Myr. (i) Same as (h) but for a source activity duration of 100 Myr. 



is that for injection of a proton at p = po « Pj,„t at time f=0, 
and acceleration for time t - xjc, the proton has a rigidity of 
pmax ~ [-^^/.^gainCPcu^l'^Vcm- Oncc morc, the heaviest compo- 
nents are close to the default case of Fig. [8] because of their 
higher charge (and acceleration rate), even for such a short ac- 
celeration time. 

5. Conclusions 

We have studied DSA of protons and nuclei in astrophysical 
sources in the presence of background radiation fields and there 
are some general conclusions we can draw. Provided the ambi- 
ent magnetic field is not too high (see Fig. |4] for protons) we 
find that when the acceleration mechanism is not limited by 
confinement protons usually reach higher energies than nuclei 
which are more limited by interactions. This conclusion is par- 
ticularly true for the conservative assumption that the turbu- 



lence in the vicinity of the shock is of Kolmogorov-type (i.e. 
a slow dependence of the acceleration rate with energy). If the 
turbulence brings the diffusion coefficient dependence closer to 
the Bohm regime we find that heavy nuclei can usually reach 
maximum energies as high as those of protons, even in the case 
of very high ambient photon fields. Whenever the acceleration 
mechanism is not confinement-limited we find there to be no 
trivial scaling of the maximum energies of the various species 
with their mass or atomic number. 

It is beyond the scope of the present paper to calculate the 
expected spectra at Earth for each source model, and we re- 
fer the reader to Allard et al. (l2005H2008"b for detailed stud- 
ies of the propagation of cosmic-ray nuclei in the extragalactic 
radiation fields. However, the non-trivial scaUng of the maxi- 
mum energy with the mass or atomic number should not af- 
fect significantly the results obtained in propagation studies in 
which maximum energies are usually assumed to be propor- 
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Fig. 11. As Fig.|6]except that p^ax = lO'"* V. 

tional to the charge. Indeed, we have seen that secondary nu- 
cleons produced in interactions within the accelerator continue 
to be accelerated (except for escaping neutrons) and can then 
reach higher energies. On the other hand, nuclei at energies 
above ~ Ax (4x10'^ eV) are expected to be almost fuUy photo- 
disintegrated by CMB photons during propagation (see Allard 
et al.. l2008l l. and this means that nuclei accelerated above these 
energies would not be seen at Earth and would simply release 
their nucleons at energies a factor A lower Whether or not nu- 
clei above ~ A x 4x 10'^ eV are photo-disintegrated within the 
source or during the propagation will make little difference to 
what is observed at Earth. Furthermore, any bumps in the spec- 
trum due to secondary nucleons should be smoothed by source 
to source variability as well as by cosmological effects during 
the propagation. Bumps in the cosmic ray spectrum would re- 
main after propagation of cosmic -rays from a single source, 
but not remain after summing the contributions of all sources 
in the Universe because of the large range of propagation times 
from the sources to the earth except for the case of a dominant 
contribution by a very nearby source. 



One important point to emphasize is that even in the pres- 
ence of strong photon backgrounds, nuclei are accelerated to 
energies above 10''^ eV whenever protons are, even though nu- 
clei are more limited by interactions than nucleons. Hence, our 
results show that DSA at sources with physical conditions that 
would cause nuclei to be fully photodisintegrated during ac- 
celeration but allow protons to reach energies above 10^*^ eV 
are excluded. In other words a pure proton composition is un- 
likely to result from an acceleration process whenever a mixed 
composition is injected at a shock. A 100% proton composition 
would therefore require a different explanation such as injec- 
tion of only protons at the shock, or propagation in a medium 
optically thick to photodisintegration of nuclei in the vicinity 
of the source. 

It could be argued that our modeling of the photon back- 
ground in which we scaled the extragalactic IR/Opt/UV by fac- 
tors between 1 and 10000 is too simplistic. However, we shall 
now consider our results briefly in the context of acceleration 
in the radiation field at a shock at a hot-spot of an FR-II radio 
galaxy. FR-II radio galaxies are thought to be the UN-beamed 
counterparts of FSRQ. Given the spectral energy distribution 
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(SED) of a "typical" FSRQ, we may infer the radiation field 
at the likely diffusive shock acceleration sites at the far end of 
an FR-II jet. In Fig. [T] we showed the radiation field expected 
at 200 kpc along the jet from 0208-512 based on the SED for 
this FSRQ given in fig. 5 of Tavecchio et al. ( 12002b . and com- 
pared it with the IR/Opt/UV at z=0, and IR/Opt/UV fields that 
are 100 and 1000 times higher. We see that while the low en- 
ergy part of the SED is important, the high energy part of the 
SED, i.e. X-rays and above, will have little effect as target pho- 
tons for pion photoproduction or photo-disintegration of nuclei. 
However, there will actually be a very large range in the nor- 
malizations of the FSRQ blazar field at the hot-spot for two rea- 
sons. First, the distance from the core to the hot-spots for FR-11 
is probably in the range 50-800 kpc, giving rise to a range of a 
factor of about 300 in radiation density. Secondly, there is also 
a large range in luminosity of FSRQ - Zhou et al. ( I2007I I found 
that log[Lv(l-5 GHz) / W Hz '] to range from 26.4 to 29.0 with 
a mean value 27.7 for 45 FSRQ. This is to be compared with 
that of 0208-512 for which log[Lv(l-5 GHz)] ^ 28.5 W Hz '. 
Nevertheless, as the shape of the SED of the extragalactic of 
photon background is not very different from the FSRQ SED 
shown in Fig. [T] our results can then give some hint of the 
influence of the photon backgrounds in this type of source. 
However, we plan to consider more specifically the accelera- 
tion of cosmic -rays in FR-11 galaxies in a future paper 

Interaction-limited acceleration processes are also expected 
to result in production of gamma-rays and neutrinos. As shown 
by Protheroe (2004), the position of the neutrino peak can help 
to constrain the acceleration parameters in proton acceleration 
sources and the characteristics of the source. However, as can 
be seen in the various figures in Sect. 4, when we include nu- 
clei neutrino fluxes cannot give any strong constraints on the 
injected composition at the shocks. 

When the size of the source becomes comparable to or 
smaller than the Larmor radius of protons at the highest en- 
ergies, the acceleration becomes limited by confinement, in 
which case heavy nuclei are accelerated to higher energies than 
protons unless the acceleration is very slow and/or the pho- 
ton background is very high. In the case of the maximum pro- 
ton energy being below 10^" eV the maximum energies of all 
species are proportional to their charge, even in the presence 
of high photon backgrounds, provided the acceleration param- 
eters are favorable (p'^J > 10^' eV and 6 close to "Bohm scal- 
ing"). Such a composition is especially interesting in the con- 
text of recent cosmic-ray data. Indeed, preliminary results from 
the Pierre Auger Collaboration show the composition possi- 
bly getting heavier above 10''' eV (see linger et al., 2007). If 
this trend is confirmed, the most likely explanation would in- 
volve confinement-limited acceleration with the maximum en- 
ergy of the species scaling with the charge. A composition get- 
ting heavier above lO'^ eV is indeed very difficult to obtain 
with propagation effects only (see Allard et al., 12005112008b 
when protons are accelerated up to the highest energies. We 
have also seen in previous sections that in most cases if the 
acceleration is interaction-limited then nuclei are not acceler- 
ated to higher energies than protons. Interestingly, it has been 
shown recently (Allard et al., ,2008, ) that the highest energy 
cosmic-ray spectrum can be successfully fitted by assuming 



a low proton maximum energy (~ 10'^ eV) and a charge- 
scaling maximum energy for the other species, providing abun- 
dances of heavy nuclei slightly higher than in the Galactic cos- 
mic rays. Although other possibilities are not ruled out, such 
a scenario where protons are not required to be accelerated 
above 10^" eV but where the highest energy cosmic-rays are 
provided by the heavy component would completely change 
the expectations and constraints on possible accelerators like 
blazars (Miicke et al., 2003) or radio-galaxies. We shall in- 
vestigate these implications in forthcoming papers. Finally, we 
note that confinement-limited acceleration is unlikely to pro- 
vide strong neutrino fluxes at high energy. 
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